Mixing-induced CP violating sources for electroweak 
baryogenesis from a semiclassical approach 
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Abstract. 

The effects of flavor mixing in electroweak baryogenesis is investigated in a generalized semiclassical WKB approach. 
Through calculating the nonadiabatic corrections to the particle currents it is shown that extra CP violation sources arise from 
the off-diagonal part of the equation of motion of particles moving inside the bubble wall. The mixing-induced source is of 
the first order in derivative expansion of the Higgs condensate, but is oscillation suppressed. The numerical importance of 
the mixing-induced source is discussed in the Minimal Supersymmetric Standard Model and compared with the source term 
induced by semiclassical force. 
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Electroweak baryogenesis (EWBG) is a promising 
scenario for explaining baryon number asymmetry in 
the universe, unlike other scenarios valid at grand uni- 
fication scale it can be tested by the upcoming collider 
experiments. An efficient way to generate large baryon 
number asymmetry is to generate it nonlocally, through 
a charge transportation mechanism [1]. In some mod- 
els such as MSSM and the two-Higgs-doublet model, 
the bubble wall of the electroweak phase transition is 
typically thick 13, ^ for the particles giving dominant 
contributions, namely the particle typical Compton wave 
length A ~ 1/r is much shorter than the wall width 
Ln- ^ (10 — 20)/r where T is the critical temperature. 
In this regime the validity of semiclassical approach 
M Isl la LZI ISD should be justified for single flavor case, 
which provides an intuitively simple description by treat- 
ing particle transportation as motion of WKB wave pack- 
ages. In a slowly moving and CP violating Higgs conden- 
sate background, the dispersion relation for particles and 
antiparticles are modified differently, contributing to dif- 
ferent semiclassical forces, which leads to a net excess 
or deficit of the particle number which can be converted 
into left-handed fermion number asymmetry. The asym- 
metry of the local fermion density then get transported in 
front of the bubble wall, which bias the baryon number 
violating processes. 

The problem becomes more involved in the multiple 
flavor case, where the CP violating mass matrix has non- 
trivial space-time dependence. The CP violating effects 
may show up not only in the dispersion relations, but 
also in the mixing of states. In this work we general- 
ize the semiclassical method by taking into account the 
non-adiabatic corrections from the spatially varying fla- 
vor mixings [9]. We restrict ourselves in the parameter 



region where the mixing between the local mass eigen- 
states are relatively small, and can be treated as perturba- 
tions to the equation of motion for local mass eigenstates. 
In two flavor mixing case the explicit form of the equa- 
tion array for left-handed component can be rewritten as 
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Since the off-diagonal elements Di2{2i) contain differen- 
tial operators, one can not obtain decoupled equations for 
4i(2) separately. 

From the dispersion relation one can deduce the group 
velocity Vg = {dco/dpc)z and semiclassical force Ft = 
(o{dvg/dt) respectively 
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where co is the energy of the state, E — U^dU and A = 
U^MdM^^dU with M the mass matrix in flavor basis 
and U is the rotation matrix diagonalizing M^M. pc is 
the canonical momentum and Pq = (0^— nt^. Note that 
only the second term in the force term is CP violating, 
and is proportional to the spin of the local mass states. 



Taking the off-diagonal terms D12 and D21 as pertur- 
bations, the solutions can be written in a generic form 
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where D^j is the lowest order solution satisfying 
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and l\' are the corrections due to the off-diagonal terms. 

The lowest order solution L]^- is obtained by the usual 
WKB wave ansatz 

4o) = Vf/'"(^'M.', 



FIGURE 1. (a)Left, the semiclassical force induced source 
term as a function of z (GeV^ ' )■ (b) Right, the mixing-induced 
source term as a function of z- The two curves in each plot cor- 
responds toL,v — lO/T (solid) and 15 /T (dashed) respectively. 
The MSSM parameters are fixed at M2 = 150 GeV and /X = 100 
GeV, with <j)^ = 0.02. 



where pd is the canonical momentum and the function 
w, provides the correct normalization for L^,- . 

Substituting the off-diagonal terms in to the equation, 
the first order perturbation takes the following form 
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which are mixtures of the two unperturbed states. The 
mixing parameter for particle / to the first order of deriva- 
tive is given by 
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The mixing parameter for particle j can be simply 
obtained by replacing / <-> j from the above expres- 
sions. It is clear that the expansion is only valid for 
d^nii / [Anp-) <C 1. The momentum dependencies comes 
from the differentiation operators in off-diagonal ele- 
ment Di2. 

The CP violating force term, as it is proportional to 
spin s, only contribute to the spin-weighted density. To 
facilitate the comparison with the force term, we give the 
spin-weighted mixing-induced source term 
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which is a momentum odd function. 

Following the standard process to derive the trans- 
portation equations, one arrives at the usual form of the 
diffusion equation 
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where IT is a similar quantity as E but for right-handed 
field. The functions g^ ^ are given by 
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with Tfi^ 



- ^,r^, the interaction rate. 



The mixing and force induced source terms have the 
following form after thermal-average 
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where ^ = /i,/r the rescaled chemical potential and 
D, = {{p^/co)'^) /T'j the diffusion constant. 

In MSSM, the chargino transportation provides a dom- 
inant CP violating source to electroweak baryogene- 
sis. The asymmetry in chargino number is converted 
into the asymmetry in left-handed top quarks through 
Yukawa interactions. The two Higgsino SU{2) doublets 
are hi = {h^^,h^j^)^ and J12 — (hii^^^Lf respectively. 
Together with the two charged gauginos W^ and W[ , 
the charginos are combined into two left- and right- 
handed four component spinors as i/^ = {Wj^ ,h2if and 
Wr = {{^lY ^ {^Hil'^y ■ '^^^ chargino mass term has the 
form \jf]fM{z)nfL in the wino-higgsino space with mass 
matrix 
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where M2 and jj. are soft supersymmetry breaking pa- 
rameters containing CP phases and Hi (z) and H2{z) the 
Higgs vacuum expectation values (VEVs). 

In FiglT] we give the source term calculated from 
Eq.(|9]l and rescaled by 1 /T from both semiclassical force 
and mixing-induced source term for wall width L„. — 
10/r and 15/r respectively. The MSSM softbreaking 




FIGURE 2. The same as Fig[I]with M2 = 250 GeV. 



parameters are fixed atM2 = 200GeV and |^| = lOOGeV 
corresponding to a chargino mass differences m^ — mj = 
111 GeV in the broken phase. The CP phase 0^ is set to 
a typical value of (^^ = 0.02. The curves given in FiglT] 
show that both of source terms have nontrivial spatial 
dependencies. However, their origin are quite different. 
The variation of the force term comes from the third 
derivative of the kink-type Higgs condensate, which has 
one minimal and two maximums. While the mixing- 
induced source term varies due to both the wall profile 
variation and the oscillation. The oscillation leads to 
multiple local minimum appearing in the curve and is 
suppressed at large distance by the wall profile. For 
M2 around 200 GeV, the mixing induced source term 
peaks at z ~ 0.03 with an amplitude Sm/T ~ 1 .6 x 10^^, 
much larger than that from semiclassical force which 
peaks at z ~ 0.08 with Sf/T = 0.6 x 10"^. In Fig|2l 
we give the results for a larger M2 = 250 GeV, and still 
fix \lJ.\ — lOOGeV, corresponding to a larger chargino 
mass differences mi — OTj = 157 GeV in the broken 
phase. One sees that for a larger chargino mass difference 
both the source terms becomes smaller as they are 1/A 
suppressed. The oscillation in the mixing-induced source 
term becomes obvious and the wave length is shorter, 
thus the oscillation suppression is stronger The mixing 
induced source term peaks at z ~ 0.03 with an amplitude 
Sm/T ~ 8 X 10^^. Although significantly reduced, it still 
much larger than that from semiclassical force which 
peaks at z ~ 0.08 with Sf/T ~ 4 x 10^1 

To estimate the oscillation suppression effects it is 
useful to define an averaged source over the wall width 
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We calculate the averaged source for different M2 = 
200 - 500 GeV and Ust the results in Tab.l in Ref.l9J 
The results show that Sm dominates over Sp in the range 
200GeV < M2 < 350GeV. With the value of M2 in- 
creasing, the averaged source term Sm drops rapidly. For 
Ln = 10/r, atM2 = 350 GeV, the mixing induced source 
is only 5.4% of that at M2 = 200 GeV, the suppression 
is due to the increased oscillation frequency. The sup- 
pression in force term Sf is mainly from the factor 1/A, 



which makes it decrease slowly. At M2 = 350 GeV, it 
is still about 30% as large as that at M2 = 200 GeV. At 
M2 = 200 GeV, their relative size is Sm/Sf = 14.6. For 
large M2 — 350 GeV, although they are close in size, the 
mixing source still dominates with Sm/Sf ~ 2.63. This 
dominance has a mild dependence on the wall width. 
For Lw = 15/r, the relative size between the two kind 
of sources remains roughly the same, although both of 
the source term becomes smaller The mixing-induced 
source (semiclassical force) term is ~ 40% (~ 30%) of 
that at Lh, = 10/r. When the value of M2 is around 450 
GeV, the two type of source term becomes comparable 
in size. For a very large M2 — 500 GeV, the semiclassi- 
cal force term becomes dominate. The mixing-induced 
source term is about an order of magnitude smaller. 

In conclusion, we have studied the effects of fla- 
vor mixing in a generalized WKB approach in which 
the off-diagonal terms in the equation of motion are 
taken into account as perturbations. With the presence 
of a slowly moving CP violating bubble wall, an extra 
mixing-induced CP violating source appears which ex- 
hibit an oscillation behavior in analogy to the neutrino 
mixings. We have made a numerical study of the oscilla- 
tion suppression effects for chargino case in MSSM and 
shown that for a light 200 < M2 < 350 GeV even in the 
small mixing case the mixing-induced source already in- 
dicate that a significant enhancement of the final baryon 
number asymmetry is possible, which will make MSSM 
a more realistic model for EWBG. 
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